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We present the first measurement of the production cross section of a W boson with a single charm 
quark (c) in pp collisions at */s = 1.96 TeV, using soft muon tagging of c jets. In a data sample of 
~ 1.8 fb _1 , recorded with the CDF II detector at the Fermilab Tevatron, we select events with W + 1 
or 2 jets. We use the charge correlation between the W and the muon from the semileptonic decay of 
a charm hadron to extract the Wc signal. We measure owe {pre > 20 GeV/c, |7j c | < 1.5)xBR(W — > 
hi) = 9.8 ± 3.2 pb, in agreement with theoretical expectations. 

PACS numbers: 13.20.Fc, 13.38.Bc, 13.85.Lg 
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In the standard model (SM), the associated production 
of W bosons with single charm quarks (Wc) in proton- 
antiproton collisions is described at leading order (LO) 
by the scattering of a gluon with a down, strange or bot- 
tom quark. Despite the dominant d-quark parton den- 
sity function (PDF) the first process contributes only 
~ 10% at the Tevatron, suppressed by the small Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element V c d- The 
contribution from gb — » Wc is also heavily suppressed 
by V e b and the b quark PDF. Therefore, about 90% of 
the Wc signal is produced by strange quark-gluon fusion 
and the production cross section is directly sensitive to 
the gluon and s-quark PDFs in the proton [l[ , at a scale 
of order of the W mass (Mjy), and to the element of 
the CKM matrix V cs . Calculations of W + heavy quark 
production are available at LO and next-to-leading order 
(NLO) in QCD [2]. From direct measurements, \V ca \ is 
known with a precision of ~ 35% from W + — > cs, and 
of ~ 10% from D, D s decays 3], largely due to the- 
oretical uncertainties. The s-quark PDF is known from 
neutrino-nucleon deep inelastic scattering experiments at 
momentum transfers, Q 2 , about three orders of magni- 
tude lower than Myy [4j . Moreover, a lack of constraints 
on the shape of the s-quark PDF has been shown to affect 
the calculation of the Wc cross section by up to 20% @ . 
The dependence on the factorization and renormaliza- 
tion scales in the NLO calculation yields an additional 
uncertainty of ~ 20% [2|. Overall, the uncertainty on 
the theoretical expectation of Wc production is ~ 30%. 
Outside the SM, an exotic particle with high semilep- 
tonic branching ratio (BR) produced with a W boson, as 
in the hypothesis investigated in @, 0, [1| , could lead to a 
measured cross section of Wc higher than expected. In 
addition, charged Higgs production via c + s — > H + is 
sensitive the the strange PDF @. Finally, Wc is an im- 
portant component of the W + 1, 2 jet samples, a dataset 
used in the search of signals such as single top, the Higgs 
boson and supersymmetric top. In this context, the mea- 
surement of the associated production of W bosons with 
single charm at the Tevatron is an opportunity to test 
the SM calculation and set an experimental constraint to 
the theory. 

This Letter documents the first measurement of Wc 
production in pp collisions. We use the correlation be- 
tween the charge of the W and the charge of the muon 
from the semileptonic decay of the charm hadron in the 
W + 1, 2 jet sample. Charge conservation in the pro- 
cess gq — > Wc (q = d, s) allows as final states only the 
pairings W + c and W~c. It follows that the charge of the 
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lepton from the semileptonic decay of the c quark, and 
the charge of the lepton from the W decay, are always of 
opposite sign. The observed lepton charge correlation in 
Wc is expected to be less than 100%, due to dilution from 
hadronic decays in flight and hadrons misidentified as 
muons ( "mistags" ) . We identify the W by its decay to an 
isolated electron (muon) carrying large transverse energy, 
Et (momentum, px) with respect to the beamline, plus 
a neutrino. We refer to these high-py electrons or muons 
as "primary leptons" (PL) . The neutrino escapes the de- 
tector causing an imbalance of total transverse energy, 
referred to as "missing Et"($t) @- Quarks hadronize 
and are observed as jets of charged and neutral parti- 
cles. We identify the charm jets from their semileptonic 
decay by looking for a muon within the jet ( "soft-lepton 
tag" or "SLT"). We use the difference between events 
in which the PL and SLT are oppositely charged ("OS") 
and events in which they have the same charge ("SS") 
as an observable and define the "asymmetry" : A — 
(Nos~ Yss)/(iVos+-Yss)- We measure A and determine 
the expected number of OS— SS events due to background 
sources, primarily VF+ light flavor (LF) hadrons, di-muon 
events from Drell-Yan (DY), and multi-jet events (re- 
ferred to as "non-VF QCD"). To a smaller extent, single 
top, 2 — ► tt and WW events are also sources of charge- 
correlated background. Because the Wc signal produces 
mainly OS events, whereas most of the background is 
equally likely to give OS or SS events, we are able to 
extract the signal from a background-dominated sample 
using the excess of OS— SS events over the expected (not 
including Wc) number. 

The CDF II detector is described in detail elsewhere 
(To| . The components relevant to this analysis include the 
silicon vertex detector, the central outer tracker (COT), 
the central electromagnetic and hadronic calorimeters, 
the central muon detectors and the luminosity coun- 
ters. The data sample, produced in pp collisions at 
y/s = 1.96 TeV during Run II of the Fermilab Teva- 
tron, was collected between March 2002 and April 2007. 
This analysis is based on an integrated luminosity of 
1823 ± 109 pb _1 . Monte Carlo (MC) simulations of 
W + jets production are performed using ALPGEN (v2 . 1 



11 



), coupled with CTEQ5L PDFs [T3| and PYTHIA (v6.3 
_13j |) for the shower evolution. Modeling of b and c hadron 
decay is provided by EvtGen [Hj , DY and WW produc- 
tion are modeled with PYTHIA (v6.5). Events are se- 
lected with an inclusive lepton trigger requiring an elec- 
tron (muon) with E T > 18 GeV (p T > 18 GeV/c). Fur- 
ther selection requires that candidate PLs are isolated 
(I < 0.1 [H) and have E T (p T ) greater than 20 GeV 
with |^| < 1.1. The event must also have $ T > 25 GeV. 
Jets are identified using a cone algorithm [l6| with a cone 
opening of AR = (A-q) 2 + (A0) 2 = 0.4 and are con- 
strained to originate from the pp collision vertex. Muons 
inside jets are identified by matching the tracks of the jet, 
as measured in the COT, with track segments ("stubs") 
in the muon detectors. Such a muon with px > 3 GeV/c 
and within AR < 0.6 of a jet axis is an SLT. The prob- 
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ability of misidentifying a hadron as an SLT muon, de- 
noted as the "mistag probability" , is measured using jets 
from data [13] ■ An SLT mistag probability look-up ta- 
ble is constructed as a function of the track px, r\ and 
(j). The number of tags predicted with the look-up table 
agrees to within 10% of the observation in a variety of 
samples, including QCD multijets and Z + jet events. To 
reduce background from di-muon resonances and double- 
semileptonic B hadron decays we remove events in which 
the muon PL and SLT are oppositely charged and have an 
invariant mass consistent with a Z, T or, irrespectively 
of the PL flavor, less than 5 GeV/c 2 . Remaining cuts are 
as in [13] • The jet energies are corrected to account for 
variations of the detector response in rj and time, as well 
as multiple pp interactions, but are not corrected back to 
the parton energy. Finally, the W ± jet dataset is parti- 
tioned according to the number of jets with Et > 10 GeV 
and \rj\ < 2.0 in the event (corresponding to the charm's 
l^cl ^ 1-5, from the geometric acceptance of muons, and 
pre ^ 20 GeV/c). We refer to events in which the PL 
is an electron (muon) as the "electron sample" ("muon 
sample"). In the electron (muon) sample we select 636 
(425) OS and 450 (313) SS events with at least one jet 
tagged by the SLT algorithm. In rare cases where there 
is more than one tag in the same jet, or more than one 
tagged jet in the same event, we use the SLT candidate 
that has the best match between the COT track and 
muon stubs. 

The composition of the W + 1, 2 jet sample includes 
primarily W_ + LF, Wbb, Wcc, Wc, non-W QCD, and DY 
events. Wbb and Wcc backgrounds reduce the statistical 
precision of the measurement but do not contribute to the 
measured asymmetry, since the presence of both heavy 
quark and antiquark allows with equal probability either 
charge combination of the leptons from the W and from 
the (anti)quark decays. Table U shows the data and ex- 
pected background contributions. The dominant mecha- 
nism by which DY events contribute to the background 
is when one of the produced muons emits a high-energy 
photon, and the muon-photon pair is identified as a jet, 
while the muon passes the SLT-tagging requirements. In 
order to reduce this background we reject events based 
on the jet's electromagnetic energy fraction and number 
of tracks. Events with an OS primary muon and SLT 
that have a jet with less than 2 tracks in a cone of open- 
ing AR — 0.6 around the jet axis and an electromagnetic 
fraction higher than 0.8 are rejected as potential radiative 
DY events. After the rejection we estimate the number 
of remaining DY events in the sample by measuring the 
number of OS-SS events in the data in a window of in- 
variant mass containing the Z peak, n|- ata , and using MC 
simulation to extrapolate n| ata in the Wc signal region. 
The uncertainty shown is due to the size of the data and, 
to a lesser extent, MC samples. 

Events without W bosons that enter the signal sample 
are typically QCD jet events where one jet has faked a 
high-j»T lepton, mismeasured energies produce apparent 
fx an d an additional jet contains an SLT muon. A frac- 



tion of these events are from bb and cc, where the PL re- 
sults from a semileptonic decay on one of the fragmenting 
heavy quark and the SLT from a semileptonic decay on 
the other, with a potentially large charge asymmetry. We 
estimate the number of non-W QCD events in the sample 
by extrapolating the number of SLT-tagged events with 
an isolated lepton (/ < 0.1) and < 15 GeV into the 
signal region of $ T > 25 GeV. The difference seen in 
Table U between the electron and muon channels reflects 
the higher likelihood for a jet to be misidentified as an 
electron than as a muon. The extrapolation uncertainty 
is determined by comparing with data the prediction for 
the number of events in the region of 0.1 < I < 0.2 and 
$ T > 25 GeV, just outside the signal region. We assign 
a ±20% systematic uncertainty based on the difference 
between the two. We model the charge correlation of 
non-IV events from the data using events with $ T > 25 
GeV but with a non- isolated PL, and the dominant un- 
certainty is from the limited data sample size. 

W + LF events enter the data sample when one of the 
LF tracks in the jet is misidentified as a muon. Since the 
same process that describes Wc production may be ex- 
pected to describe Wu, we expect a small anti-correlation 
between the charge of the W and the charge of the tracks 
in the jets recoiling against the W. We estimate the 
total W + LF background using the mistag probability 
parametrization with the number of events before SLT 
tagging, correcting for a fraction due to non-W QCD. 
The uncertainties include contributions from the mistag 
prediction and from the correction for non-W QCD. We 
determine the expected asymmetry for these events using 
ALPGEN MC events; the asymmetry is lower in the muon 
sample as a result of the removal of events consistent with 
di-muon resonances, which removes only OS events. Re- 
maining small backgrounds with expected asymmetries 
are from single top (i-channel) production, WW, and 
Z — ► tt. These are estimated using theoretical cross sec- 
tions and MC simulations. Finally, for Wbb and Wcc we 
measure a charge asymmetry from MC events, compati- 
ble with zero as expected. 



TABLE I: Summary of data and charge-asymmetric back- 
ground in the SLT-tagged W + 1, 2 jet sample. 



Source 


Events 


Asymmetry 


OS-SS events 


Drell-Yan 






41.6 ± 5.4 


Non-W (e PL) 


200.0 ±41.5 


0.179 ±0.046 


35.8 ± 11.8 


Non-W (fx PL) 


30.2 ±6.5 


0.252 ± 0.066 


7.6 ±2.6 


W + LF (e PL) 


695.5 ± 75.7 


0.057 ±0.002 


39.6 ± 4.9 


W + LF {fj, PL) 


491.4 ± 53.3 


0.019 ±0.002 


9.3 ±2.0 


Single top 






7.6 ± 1.1 


Z->TT, WW 






7.2 ± 1.2 


Total Background 






148.7 ± 15.4 


Data 


1824 




298 



After background subtraction we measure 149.3 ± 
42.1(stat) ± 15.4(syst) Wc events in the selected SS- 
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FIG. 1: SS-subtracted distributions of the SLT muon pr rel- 
ative to (a): the beam axis and (b): the SLT-jet axis. The 
Wc contribution is normalized to this measurement. 



subtracted sample. Figure [TJ shows two kinematic dis- 
tributions of the SLT muons in the tagged events. Each 
of the distributions is SS subtracted and the stacked his- 
tograms show the contribution from both the expected 
Wc signal and the estimated backgrounds. The higher- 
Pt muons in both plots arise primarily from DY events 
misidentified as charm decay. The Wc production cross 
section is obtained using the cross section formula: 



a Wc x BR(W -> tv) = 



N, 



os-ss 



OS-SS 
bkg 



Acc ■ J L dt 



(1) 



where awe = ctw+e + &w-c with px c > 20 GeV/c 
and \rj c \ < 1.5, and BR(W -> tv) = 0.108. In the de- 
nominator, Acc is the SS-subtracted acceptance times 
efficiency derived from a MC simulation of Wc events, 
and J Ldt is the integrated luminosity of the data. We 
find Acc = (0.833 ±0.010) x 10~ 2 ; the value implicitly in- 
cludes the expected charge asymmetry of the Wc sample, 
which is (70.6 ± 0.5)%, as well as the semileptonic BR of 
charm hadrons to muon. The uncertainty is due to the 
size of the MC sample. Using Equation[TJ we find: awe X 
BR {W -> tv) = 9.8 ± 2.8(stat)t i ;g(syst) ± 0.6(lum) pb. 



TABLE II: Summary of systematic uncertainties. 



Source 


Aa Wc (%) 


Lepton ID 


±1.4 


SLT tagging efficiency 


±5.1 


Jet energy calibration 


+4.0/-2.7 


PDFs 


±3.6 


c-quark hadronization 


±4.6 


ISR, FSR 


F4.2/-9.5 


Factorization scale 


±1.3 


Background 


±10.3 


Luminosity 


±6.0 


Total (excluding Luminosity) 


+14.2/-16.3 



Systematic uncertainties are shown in Table [TTJ MC 
modeling of the efficiency for identifying the PL ("lep- 
ton ID") is measured using Z-boson data and MC sam- 
ples. The uncertainty on the efficiency for identifying the 
SLT includes a contribution from MC modeling of par- 
ticle tracking efficiency in dense jet environments, and 
a contribution from the parametrization of the SLT effi- 
ciency as a function of the px of the muon [l7J . The PDF 
uncertainty is derived by using CTEQ and MRST sets, 



as prescribed in 1S|. We compare charm jets in PYTHIA 
and HERWIG [l|| 2C | to evaluate the uncertainty due to 
different hadronization models. To measure the effects of 
initial (ISR) and final state radiation (FSR) we use inclu- 
sive Wc production in ALPGEN and compare the nominal 
acceptance with enhanced or reduced radiation as in [l8| . 
The uncertainty on the factorization and renormalization 
scales is estimated by varying them in ALPGEN between 
the transverse mass of the W and that of the charm. 
Finally, uncertainties on the background estimations are 
included in the cross section as systematic uncertainties. 

In conclusion, we present the first measurement of 
the Wc production cross section using ~ 1.8 fb -1 
of pp collisions at y^i = 1.96 TeV. We mea- 
sure a Wc (p Tc > 20 GeV/c, |t? c | < 1.5) x BR (W -> tv) = 
9.8 ±3.2 pb, in agreement with the NLO calculation over 
the same phase space region of: ll.Oj^'o pb [2l|. The 
result provides an experimental validation of the theo- 
retical Wc cross section for use in single top and Higgs 
search analyses. Within the quoted sensitivity, we see no 
evidence of the presence of an exotic particle with high 
semileptonic BR. Finally, the result provides the first ex- 
perimental constraints to the theory uncertainties; an ex- 
trapolation indicates that a precision of ~ 15% can be 
achieved with ~ 6-7 fb" 1 of data, expected by the end 
of Run II. 
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